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The frequencies of quasi-periodic oscillations around neutron stars are believed to be related to 
characteristic frequencies in the gravitational fields of the compact objects. In different variability 
models, these include the Keplerian, epicyclic, and Lense-Thirring frequencies, which depend mostly 
on the properties of the stellar spacetimes. We argue that quasi-periodic oscillations in the X-ray 
flux of neutron stars can be used to map the external spacetimes of the compact objects and, 
therefore, lead to direct tests of general relativity in the strong-field regime. In particular, we show 
that particular extensions of General Relativity, in which the gravitational force felt by matter is 
mediated by both a rank-two tensor and a scalar field, can be constrained by current observations. 



I. INTRODUCTION 

Neutron star properties, such as radius and mass, help 
probe the nature of strong-field gravity. In a previous pa- 
per pj, we showed how the measurement of the redshift 
of atomic lines emitted from the surface of a neutron star 
can be a powerful tool for investigating the properties of 
gravity in the strong-field regime. In particular, we found 
that our current lack of knowledge about the behavior of 
gravity in the strong-field regime far outweigh our cur- 
rent lack of knowledge about the neutron star equation 
of state. 

In this paper, we expand that investigation to show 
how the properties of the spacetime around a neutron 
star are at least as sensitive to the choice of gravity the- 
ory in the range currently allowed by weak-field observa- 
tions as they are to the exact nature of the neutron star 
equation of state. As in 0], we choose a particular single- 
parameter scalar-tensor theory but this time we extend 
our consideration to the spacetime outside the neutron 
star surface. 

In particular, we examine the way in which devia- 
tions from general relativity may affect the frequencies 
of quasi-periodic oscillations (QPOs) found in many X- 
ray binaries. QPOs are of contested origin; we refer the 
reader to Refs. 0,0. All models, however, agree that 
the processes responsible for QPO production take place 
in a region of order the Schwarzschild radius around the 
compact object. This makes QPOs excellent probes of 
strong-field gravity, especially if, as many believe, their 
frequencies are directly related to the characteristic fre- 
quencies of test particles near the innermost stable circu- 
lar orbit (ISCO). The growing numbers of QPO observa- 
tions promises an ideal testing ground for the predictions 
of general relativity. 

We provide two new methods of constraining devia- 
tions from general relativity with QPO observations. The 
first method relies only on general arguments about the 
maximum coherent frequency of the lowest-order modes 
near the surface of neutron stars [j] . The second method, 
which relies upon a particular interpretation Q of a 



commonly observed QPO phenomenon - the correlation 
in peak position found between high frequency QPO 
pairs 0, 13 ~ approaches the binary pulsar tests in its 
ability to constrain the particular scalar-tensor effect con- 
sidered. The generality of our first method demonstrates 
the importance of examining a wide range of high-energy 
X-ray phenomena of diverse nature and origin for clues 
about the properties of strong-field gravity. The power of 
our second method demonstrates how the development 
of more comprehensive physical accounts of these phe- 
nomena can lead to new advances in our ability to use 
the universe as a gravity laboratory, probing spacetime 
curvatures many orders of magnitude higher than those 
previously accessible to precision gravity tests. 

Our paper is organized as follows. In section two, we 
describe the field equations and stellar structure equa- 
tions for stars in a general scalar-tensor theory and we 
solve them through a combination of analytic and com- 
putational techniques. In section three, we discuss in 
detail the previously unremarked-upon properties of test 
particle orbits around such stars and, in section four, we 
apply these results to show explicitly how QPO proper- 
ties can be used to constrain the nature of strong-field 
gravity. 

II. EQUATIONS OF STELLAR STRUCTURE 

Following |8( as an example of a gravity theory that 
deviates from GR in the strong-field regime, we consider 
a theory containing a scalar field coupled to gravity with 
the action 

S = yr^r J d^x^iR* - 2g^^A,u) 

+S m [$m,A 2 ((l>)g* tlv ]. (1) 

Here fy m refers collectively to all matter fields other than 
4>, G* is a dimensional constant, and A((f>) is a function of 
4>. The (0, 2) tensor g* M „ is the Einstein frame metric; all 
matter (apart from </>), flows on geodesies of the physical 
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(or "Brans-Dickie") frame metric, A 2 {<jj)g* l]bV , which we 
will refer to as g M „. The field equations for the scalar field 
(f> are particularly simple in the Einstein frame where <j> 
appears to be non-minimally coupled to gravity; because 
of the unphysical nature of the metric g*^, however, 
confusion may arise; the reader is referred to Q for a 
deeper discussion of these issues. Unless otherwise spec- 
ified, all quantities reported in this paper are measured 
in the physical frame. 

We summarize the field equations, following [Icj : 



ii* M „ — 2d ll 4>d v 4> + 87rG* [t*^ — tT*q*iiv 
V»„V^ = -47rG*a(<£)T*. 
The equation of motion for matter is: 

V*!/T^ = a(</>)T»V» M 0, 

where starred subscripts refer to operations using the 
Einstein frame metric, and a{<j>) = d\nA((f>)/d(f>. The 
stress energy tensor in the Einstein frame, T* , is equal 
to A & (<f))T^ v , where T^ v is measured in the physical 
frame. 

We see that the scalar field, </>, is sourced by the trace of 
the stress-energy tensor, proportional to the logarithmic 
derivative of the warp factor, A{4>). Thus, though the 
cosmological value of <fi may be zero, a strong enough 
concentration of matter can support a non-zero value of 
4>. Given a sufficiently dense concentration of matter, and 
a sufficiently steep A(<f>), a region of space with non-zero 
<j> may become energetically favorable. 

Solar system experiments constrain the function A to 
be very flat at the cosmological value of <fi, denoted by 
4>q. If, however, the <p field inside a compact object fluctu- 
ates just far enough away from cosmological to discover 
a steep part of A, perhaps during a violent formation 
event, the system will be able to reach the more energet- 
ically favorable configuration, with large non-zero values 
for the scalar field </> near the center of the object. This 
is the phenomenon of "spontaneous scalarization." 

This phenomenon occurs in general for functions A(<p) 
when the derivative din A/difil^g is sufficiently nega- 
tive ^} 1 although the precise configuration of the neu- 
tron star and its external spacetime may vary with the 
particular choice of A{(j>). We consider only the first term 
in the Taylor expansion of In A{<j>) that leads to negative 
curvature, 



A(<p) 



e^ 2 



(2) 



where (3 is a real number. This is the same form as 
used in where it is referred to as the "quadratic 
model" because of its relationship to the Taylor expan- 
sion of In A{4>). Other investigations of scalar-tensor the- 
ory have chosen different forms for A{<j>); see for a 
comparative list of choices of A(<j>), as well as for how to 



rewrite the action in Eq.^cxplicitly for particular choices 
of A(4>) used in the literature, in both the Einstein and 
physical frame. 

Again following we can write the relativistic equa- 
tions of stellar structure for a neutron star. We write our 
metric as: 



dst 



^dt 2 



Mr) 



dr 2 



+r 2 (d9 2 + sin 2 ed(f> 2 ), 



(3) 



and describe the matter fields as a perfect fluid in the 
physical frame: 



T. 



(p + pju^Uv+pg^u. 



The differential equations to be integrated are then: 
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Here N is the baryon number, n the number density, 
p the density and p the pressure. We supplement these 
equations with a (zero temperature) equation of state, 
p = p(p) and n = n(p). Primes denote derivatives with 
respect to coordinate radius r. 

Having chosen an equation of state, we then integrate 
Eqs.QJjniirom the center of the star, where we specify 

p(0) = i/(0) = 0, 
0(0) = 0c, 
P(0) = Pc, 

N(0) = n c , 
p(0) = Pc, 

to the surface, where p = 0. We then must integrate the 
equations for is, (j), tp, and p from the surface of the star 
to infinity to determine the form of the metric exterior 
to the star. 

Throughout this paper we will be interested only in 
solutions where <f>o, the cosmological value of <fi, is zero 
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[and A((f>o ) = 1, which can be achieved by renormalizing 
A]. We create our set of solutions to Eqs. (0J-@ with 
a simple shooting method We note that \i does not 
have compact support; outside of the star the scalar field 
(j) is non-zero (though relaxing as 1/r to zero) and con- 
tributing to the physical Arnowitt-Dcscr-Misner (ADM) 
energy felt by an observer far away. In our terminology, 
the physical ADM energy is that that would be inferred, 
e.g., from test particle orbits near infinity, and is identi- 
cal to the M^DM m 13 • When considering solutions to 
the metric, we renormalize the lapse function, exp \y(r)\ 
so that the physical metric, Eq. J3J), goes to rj^, the 
Minkowski value, at infinity. 

In an earlier paper Q , we considered three commonly 
used equations of state, which cover a broad subset of the 
wide range discussed in Cook et al. Ihi. In order of in- 
creasing stiffness, these were EOS A [ll, EOS UU E3, 
and EOS L Neutron star models computed with 

these three equations of state bracket the uncertainty in- 
troduced by our inability to calculate from first principles 
the properties of ultra-dense matter, when condensates or 
unconfined u-d-s quark matter is not taken into account. 

For this paper, we will use the "median" EOS UU 
to illustrate the qualitative behavior of the variables we 
wish to examine; the results calculated for UU differ only 
slightly from those calculated with the two bracketing 
equations of state, A and L, in accord with the results 
presented in 0. 

For simplicity, we neglect the effect of rotation on the 
frequencies of orbits in the external spacetime. In Gen- 
eral Relativity, the effect of rotation on the orbital fre- 
quencies is less than 20% for the inferred spin frequencies 
of the neutron stars that show QPOs. We expect rota- 
tion to affect the conclusions in this paper at the same 
level. 



III. PROPERTIES OF SCALAR STARS 

Our gravity theory is described by the action in Eq. 
and the definition of the function A((j>) in Eq. (0). The 
function A((f>) has one free, real parameter, /?, which char- 
acterizes the strength of the scalar field coupling to grav- 
ity. 

Spontaneous scalarization occurs only when (3 < 
—4.85. In Fig. ^ we plot the binding energy, (Nrrib — 
A/adMi where m& is the baryon rest mass), as a func- 
tion of central density p c , using equation of state UU, 
for a range of values for the parameter (3. We draw the 
reader's attention to the very high binding energies that 
can be achieved when (3 is sufficiently large. This large 
energy is due mostly to energy released when <f> drops 
away from its cosmological value; when a scalarized star 
is produced, this energy goes into producing a monopolc 
"gravitational wave" in <f>. 

After the "turnover point", i.e., the central density for 
which the mass for a particular choice of (3 is at its max- 
imum, the stars are unstable to radial perturbation and 
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FIG. 1: Binding energies of our neutron star configurations, 
as a function of central density, for different values of the 
parameter j3, for equation of state UU. Dashed lines indicate 
where a star is unstable to radial perturbations. 
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FIG. 2: The mass-radius relation for our neutron star config- 
urations, for different values of the parameter 0, for equation 
of state UU. Dashed lines indicate where a star is unstable 
to radial perturbations. A fiducial value comparable to the 
most accurately measured neutron-star masses is drawn as a 
horizontal heavy dashed line. 



will collapse; these stars are marked with a dashed line in 
Fig. n Following the stability analysis of , we know 
also that when the binding energy for a scalarized star is 
greater than that for its unscalarized, General Relativis- 
tic counterpart with the same central density, only the 
scalarized configuration is stable. Scalarization is manda- 
tory when the scalarized star has greater binding energy. 

Fig. |21 shows the mass-radius relationship for the same 
range of the parameter (3 and the same equation of state, 
UU, as above. Immediately apparent is the result 0, 
that scalarized stars can achieve large masses and large 
radii. 
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IV. ORBITS IN THE VICINITY OF SCALAR 
STARS 



2.0 



Here we consider the properties of a test particle in a 
circular orbit outside the star. We first consider the or- 
bital frequency of the particle as seen at infinity. Paral- 
leling the analysis for circular orbits in the Schwarzschild 
metric [l9(, we find: 



.5 - 



1 .0 




drV (1-2/i/r) 
Ifr) ~ A 2 



E 2 



A((j)) 2 e v 



- 1- 



A(<py 



where E and / are the energy and angular momentum 
per unit mass of the test particle, respectively, r is the 
proper time of the particle, and r is the orbital radius. 
We note that, in contrast to the General Relativistic case, 
fi is a function of radius (as is, necessarily, A, which is a 
function of d>). For the orbit to be circular, we must have 
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FIG. 3: The ratio of orbital to "Keplerian" frequencies for a 
star with Mj^j)y[ = 1.5Mq for a range of values of /3, using 
equation of state UU. Not all the orbits for which frequencies 
are plotted here are stable; see Fig. [I] 
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Using these conditions to eliminate E and and using 
the lapse function dt/dr, we find the orbital frequency of 
a test particle, as measured at infinity, to be 




2(3<jxt>' 



+ 2r 



e v {l - 2fi/r)~ 1 + tp(/3(f> + rip / 2)r 2 e u / ^ 



1 + ipr(3<fr 



where, as before, primes denote derivative with respect 
to r. The final form of this equation demonstrates how 
the orbital frequency approaches the general relativistic 
value as the scalar coupling, is taken to zero. In GR, 
ip = </) = and e v (l — 2/i/r) -1 = 1, so that we obtain 
the Kepler frequency u> 2 = [i/r 3 . 

Figure [31 shows the ratio of orbital frequency to the 
"Keplerian" value as a function of the logarithm of 
the distance from the center of mass, for a star with 
A/ADM = 1.5M Q , for equation of state UU and for a 
range of values of the parameter (3. We define the Kep- 
lerian frequency at radius r to be the orbital frequency 
that would be produced at radius r were the scalarized 
star replaced by a General Relativistic star with Mj^y)M 
equal to fJ.(r). We remind the reader that (i(r) is increas- 
ing even outside the star, since </> is non-zero near the 
surface. 

Comparison of figures |3| and [3 reveals an interest- 
ing and counter-intuitive effect of spontaneous scalariza- 
tion. Scalarized stars have radii much greater than their 



General Relativistic counterparts with the same A^ADM- 
One would thus presume that their surface gravity, and 
thus their orbital frequencies near the surface, would be 
lower than the General Relativistic case. In fact, the case 
is quite the opposite; near the surface of a scalarized star 
the orbital frequency can be up to 50% higher than they 
would be at the same distance for a General Relativistic 
star with the same enclosed mass. At distances greater 
than ss 300 Schwarzschild radii, however, the orbits con- 
verge rapidly to their General Relativistic values. This 
is a dramatic illustration of how large strong-field effects 
can effectively hide themselves in a small region. 

To investigate the stability of these orbits to radial 
perturbations we calculate the epicyclic frequency k. We 
find 



1 (dr\ d 2 (dr 



2 \di) dr 2 \dr 



The orbital stability criterion is simply that k 2 > 0. The 
form of the right hand side of this equation is compli- 
cated, but it depends only on the variables in Eqs. Q) © 
and their derivatives. 

Figure 0] shows the behavior of the epicyclic frequency 
as a function of the logarithm of radius, again with 
MaDM = I.5M Q , equation of state UU, and a range 
of values for the parameter (3. A negative epicyclic fre- 
quency is used here to indicate an imaginary n (i.e., 
k 2 < 0, instability.) We see again how the properties of 
the scalarized stars are altered by factors of order unity 
compared to the General Relativistic case. The point of 
instability is changed such that, as the parameter (3 is 
made more negative, the ISCO radius is pushed closer 
to the surface of the star compared to the GR ISCO ra- 
dius. As expected, for radii much larger than the ISCO 
or stellar radius, the ratio of epicyclic to orbital frequency 
converges to the General Relativistic value. 

As can be seen in Fig. 0] for values of the parameter 
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FIG. 4: The ratio of epicyclic to orbital frequencies for a range 
of values of [3 and Af,\DM = 1.5M0, using equation of state 
UU. A negative epicyclic frequency is used here to indicate 
instability (i.e., k 2 < 0.) The limiting GR case is plotted as 
a dashed line. 

(3 < — 8 and neutron star masses -^adm ~ 1.5M©, all 
orbits down to the stellar surface are stable. 



V. CONSTRAINTS ON [3 FROM 
QUASI-PERIODIC OSCILLATIONS 

The predicted differences in orbital and epicyclic fre- 
quencies between the General Relativistic and scalar- 
tensor theories discussed in the previous section suggest 
a number of possible tests that can be conducted using 
the observed properties of quasi-periodic oscillations in 
X-ray binaries. In this section, we will discuss two such 
methods. The first, relying only on the difference be- 
tween the General Relativistic and scalar-tensor orbital 
frequencies, will turn out to be an as-yet inconclusive - 
though rather general - test. The second, which relies on 
a proposed theoretical picture of QPOs, generates tighter 
constraints. 



A. Limiting j3 by orbital frequency alone 

In many QPO models, a strong peak in the Fourier 
spectrum of the observed luminosity time-series is pro- 
duced near the orbital frequency of the inner edge of 
the Keplerian disk 3]. In both General Relativistic 
and scalar-tensor theories, the orbital frequency increases 
monotonically with decreasing distance from the center 
of the object. 

Thus, assuming that the QPO frequencies come from 
some radius close to but never less than that of the inner 
edge of the Keplerian flow where the epicyclic frequency 
becomes imaginary, this suggests a simple null test of 
scalar tensor theories. To wit: is the highest possible 
orbital frequency one can produce around a scalar-tensor 
star - for a particular value of the parameter p - lower 



FIG. 5: The highest orbital frequency produced in scalar- 
tensor theories, as a function of mass, for equation of state 
UU. The long dashed line shows the current 1330 Hz lower 
limit on the maximum QPO frequency, the short dashed line 
represents the General Relativistic case, and the dotted lines 
indicate unstable branches. 

than any observed QPO frequency? We may refine this 
question further by specifying a reasonable (ADM) mass 
range we wish the scalar-tensor stars to have. 

The current lower limit on the maximum QPO fre- 
quency is 1330 Hz [20|. We draw this lower limit on 
Fig. [SI where we show the highest orbital frequency for 
a variety of scalar-tensor theories as a function of mass, 
for equation of state UU. The nonintuitive property of 
scalar-tensor stars can immediately be seen: making the 
parameter (3 more negative increases the star radius (see 
Fig. |2J) but also increases the orbital frequency at every 
radius beyond the stellar surface (see Fig.|3|); this means 
that a wider range of masses are able to produce the 
high maximum frequency in scalar-tensor theory. There 
appears to be no simple way of demonstrating the neces- 
sity of this result; all that can be said is that the metric 
around a scalar-tensor star behaves very differently from 
its Schwarzschild cousin. 

For completeness, we present this result in a comple- 
mentary manner in Fig. [5] Here we show, as a function 
of the value of (3, the mass range for a variety of equa- 
tions of state that can reproduce the QPO lower limit 
of 1330 Hz. Thus this test, though ideal because of its 
model independent nature, cannot rule out any part of 
the parameter space. 

B. Limiting j3 by joint analysis of orbital and 
epicyclic frequencies 

The high-frequency QPOs discussed in this paper of- 
ten come in pairs |2( with the two peaks maintaining 
slowly changing separations. The relativistic precession 
model 5] predicts a simple relationship between, on the 
one hand, the positions and separations of these QPO 
peaks and, on the other hand, the orbital and epicyclic 
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FIG. 6: Allowed range (indicated by crosshatching) for 
-^ADM' as a function of j3, for EOS UU, for the current 
best lower limit on the maximum QPO of 1330 Hz. 



frequencies at the inner edge of the accretion disk of the 
neutron star. 

In this model, as the inner radius of the disk changes 
over time, the QPO frequencies change to reflect the 
change in orbital and epicyclic frequency. This relation- 
ship allows one to use measurements from a single ob- 
ject to sample, not just a single radius in the neutron 
star external spacetime, but a whole region of the orbital 
plane. This results in much stronger constraints on the 
nature of strong-field gravity while at the same time re- 
ducing the sensitivity to the choice of equation of state for 
ultra-dense matter. As we will see, this allows us, along 
with the assumption that the stars are slowly-rotating, 
to rule out the scalar-tensor theory of gravity contained 
in Eqs. Q and @ for a range of the parameter (3. 

To derive constraints on the scalar-tensor theory we 
first solve the equations of motion for particles in the 
external spacetime, for each value of the parameter (3 and 
^ADMj ^0 find the Keplerian orbital frequency, lo , and 
the epicyclic frequency, k, as a function of radius. The 
lower bound for the radius is either that of the innermost 
stable orbit, or that at the surface of the neutron star 
itself, whichever is larger. We then find the relationship 
between ui and lu — k themselves, which according to the 
model correspond to the frequencies of the two QPOs. 

We then fit this relationship to the observed QPO fre- 
quency pairs, and derive a reduced chi-squared value for 
each point in (-^ADM"^) s P ace - In Fig. [3 we plot the 
two-sigma confidence contour around the minimum x 2 
value obtained for the GR case. As can be seen, the 
onset of scalarization alters the range of allowed masses, 
and for values of (3 < —10.25, scalarization produces neu- 
tron stars with spacetimes sufficiently different from the 
GR case as to be ruled out altogether. 

We emphasize that this test has assumed that the rota- 
tion of the neutron star itself has a small (less than order 
unity) effect on the external spacetime. While we have 
not constructed and solved the equations for the scalar- 



FIG. 7: Two-sigma confidence contours (dashed line), as a 

function of -^ADM' ^ or a n ^ *° ^ ne d a t a °f R e f- EI- f° r a 
range of values of (3 and EOS UU. The regions bounded by 
solid lines show the region of standard GR behaviour ("GR" 
region) and the parameter space in which stars undergo spon- 
taneous scalarization ("S" region.) 



tensor theory in the case of a rotating matter source, we 
believe that the corrections to the orbital-epicyclic fre- 
quency relationship which governs this test will be on 
the order of 20%. 

We expect that combining observations from different 
systems, and requiring that (3 be held fixed across all 
observations, will provide additional confirmation of this 
result. Because of the uncertainty introduced by the con- 
sideration of nonzero stellar angular momentum, we can- 
not put a precise lower bound on the value of the param- 
eter (3, but given the magnitude of corrections found for 
rotation in the GR case, the lower bound is presumably 
not far from f3 ~ —10. 

The analysis of this section has relied on a particular 
model of QPO production. Other QPO models, since 
they also invoke processes at the high spacetime curva- 
tures near the surface of the star, will lead to similar 
kinds of constraints, of varying severity, for scalar-tensor 
models. 



VI. DISCUSSION 

We have investigated two different methods for con- 
straining the nature of strong-field gravity using obser- 
vations of QPOs. As an example of a theory of gravity 
that deviates from GR, we have used a particular scalar- 
tensor theory, parametrized by a single real number, (3. 
The class of theories we have investigated has gained a 
great deal of attention because of its invisibility to tradi- 
tional, weak-field tests. 

We have investigated the nature of test particle or- 
bits in the external spacetime of the neutron star, and, 
while uncovering some counterintuitive behavior, have 
confirmed for test particle orbits what has already been 
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known for other observables of compact objects in this 
scalar-tensor theory: while large (order unity or greater) 
deviations from General Relativity are observed within 
the first few hundred Schwarzschild radii of the neutron 
star, the effects of introducing a scalar field rapidly dis- 
appear at the larger distances (and far lower curvatures) 
probed by ordinary weak-field tests and require observa- 
tions of much greater precision to detect. 

The first method we have investigated, while not rely- 
ing on a particular theory of QPO production, does not 
produce strong constraints on the value of the parameter 
(3. The second method, which relies upon the relativistic 
precession model of QPO production, is more powerful 
in its ability to separate the predictions of General Rel- 
ativity from those of scalar-tensor theory. The power of 
this second method derives in part from the fact that we 
are able to sample a range of radii around the neutron 
star, and not just a single radius. 

In a previous paper [l| , we investigated the use of grav- 
itationally redshifted atomic lines in testing strong-field 
gravity, and showed that the nature of gravity in the 
strong-field is less well known than the equation of state 
for ultra-dense matter. Further constraining this equa- 



tion of state has often been taken to be the goal of theo- 
retical studies of neutron star properties. 

We find, however, that a second, promising avenue 
is opening up for investigations that attempt to con- 
strain the nature of gravity itself. In addition to the 
tests presented here and in previous work, additional 
ways of testing strong-field gravity can be performed that 
take advantage of other aspects of the rapid observa- 
tional and theoretical progress that has been made in 
recent years in the study of compact objects in the X-ray 
and soft-gamma ray bands. Another such test based on 
the Eddington-limited luminosities of radius-expansion 
bursts will be reported |22(. 

While solar system and binary pulsar tests observe 
with a much higher degree of precision than can be ex- 
pected from the investigation of high-energy processes 
near the surface of compact objects, the large spacetime 
curvatures found in the latter systems can, in some cir- 
cumstances, magnify the expected deviations and pro- 
duce constraints on the nature of GR of similar severity. 
In such cases, weak field tests and strong-field tests de- 
scribed in this paper and others may be considered com- 
plementary. 
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